Observation of Andreev bound states at spin-active interfaces 
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We report on high-resolution differential conductance experiments on nanoscale superconduc- 
tor/ferromagnet tunnel junctions with ultra-thin oxide tunnel barriers. We observe subgap con- 
ductance features which are symmetric with respect to bias, and shift according to the Zeeman 
energy with an applied magnetic field. These features can be explained by resonant transport via 
Andreev bound states induced by spin-active scattering at the interface. From the energy and the 
Zeeman shift of the bound states, both the magnitude and sign of the spin-dependent interfacial 
phase shifts between spin-up and spin-down electrons can be determined. These results contribute 
to the microscopic insight into the triplet proximity effect at spin-active interfaces. 
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In superconductors electrons are bound in Cooper 
pairs, usually in a singlet state, i.e., with opposite spin. 
Recently, flow of a supercurrent through fully spin- 
polarized chromium dioxide between two singlet super- 
conductors has been reported [l[ indicating Cooper pairs 
of equal spin, thus corresponding to a long-range triplet 
proximity effect 2]. This phenomenon is believed to be 
due to a combination of spin-dependent interfacial phase 
shifts Q of the electron wave-functions in conjunction 
with spin- flip scattering To date, most experimental 
evidence for the triplet proximity effect and spin-active 
scattering is based on the observation of the thickness 
dependence of the Josephson current [5], but there is an 
increasing interest in obtaining direct spectroscopic evi- 
dence @. In this Letter, we report on the observation of 
Andreev bound states at spin-active interfaces, allowing 
a precise determination of spin-dependent phase shifts. 

Andreev scattering at interfaces is a powerful probe 
of the superconducting order-parameter symmetry 
Zero-energy Andreev bound states have been observed 
in both singlet d-wave 18] and triplet p-wave [9] super- 
conductors. Microscopically, these bound states occur 
due to wave-mechanical phase shifts between the electron 
and the hole involved in Andreev reflection. In the case 
of unconventional bulk superconductors, a phase shift of 
TT may be induced by sign reversals of the bulk order pa- 
rameter, yielding a zero-energy bound state. Since the 
electron and the hole involved in Andreev reflection re- 
side in opposite spin bands, spin-active interfaces also in- 
troduce a phase shift 9s into the Andreev process, but in 
this case any value between —n and tt can occur. Con- 
sequently, the surface bound states can have arbitrary 
subgap energy Andreev bound states at spin-active 
interfaces are predicted to induce characteristic features 
in the den sity of states 13, Hi, in the local 13" 13 and 
non-local [13| conductance as well as in the noise spec- 
trum fl^. 

In particular, the bound states give rise to a spin- 




FIG. 1. (Color online) (a) Scanning electron microscopy im- 
age of one of our samples, together with the measurement 
scheme, (b) High-resolution transmission electron microscopy 
cross section of a reference structure, showing the atomic 
structure of the tunnel barrier. 



dependent subgap density of states, with resonances at 
the characteristic energy 



e± = ±sign(6l,)Acos(6l,/2), 



(1) 



for spin up (e+) and down (e_) quasiparticles, where A 
is the pair potential. In the differential conductance, 
the bound states show up as double-peak features due 
to resonant Andreev reflection at voltages eV ~ — e±. 
The very simple relation between bound-state energy and 
phase shift allows an unambiguous and precise determi- 
nation of Os- For small 9s, however, the features are 
expected to be very close to A, and therefore difficult 
to distinguish from the ordinary gap anomaly. Con- 
sequently, interfaces with large 9s are required for a 
clear observation of bound states. Also, the resonances 
broaden into a subgap continuum at highly transpar- 
ent interfaces, so that low-transparency tunnel junctions 
are required to clearly observe them. Recent theoret- 
ical predictions show that 9s can be of the order of tt 
in superconductor-ferromagnet junctions with ultra-thin 
tunnel barriers [l2[, and this has motivated us to fabri- 
cate this type of junctions. 
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FIG. 2. (Color online) Differential conductance spectra, (a) 
Differential conductance of one junction as a function of ap- 
plied magnetic field B at lowest temperature. The solid line 
is a fit to (13}. (b) subgap differential conductance of the same 
contact at B = (symbols), together with a fit according to 
Q (line). 



Samples with multiple junctions were fabricated by 
standard e-beam lithography and a sequence of oblique- 
angle (shadow) evaporation techniques. The junctions 
consist of a thin (« 12 — 15 nm) layer of aluminum, which 
was oxidized in situ to form a tunnel barrier, and an iron 
counterelectrode. Typical parameters of the aluminum 
films are: resistivity p ~ 10 fiflcm, critical temperature 
Tc w 1.5 K and critical field _Bc ~ 2 T for in-plane mag- 
netic fields. A total of 8 samples were prepared, with 
2 to 5 junctions each. Normal-state junction conduc- 
tances were around Gn ~ 1 mS, with a typical junction 
area of 150 x 150 nm^, as shown in Fig. [IJa). This cor- 
responds to resistance-area products as low as 30 rtp,m^ , 
and an average transmission probability of {t) w 5 x 10^^. 
In Fig. [IJb), the structure of the tunnel barrier can be 
resolved in a high-resolution transmission electron mi- 
croscopy (HRTEM) image of a reference structure fabri- 
cated by the same methods. As can be seen, the oxide 
barrier is amorphous, and varies in thickness between 1 
and 2 nm. Assuming an effective tunnel barrier height 
of about 1 eV, the tunneling probability for electrons 
through a barrier of 1 nm thickness is about 10~^, con- 
sistent with the observed resistance-area product. The 
differential conductance of these junctions was measured 
at low temperatures using a low-frequency ac method 
with an excitation of « 5 /iV at 138 Hz. Details of the ex- 
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perimental setup have been described elsewhere 

Figure [SJa) shows the differential conductance of one 
junction at lowest temperature (T — 50 mK) for differ- 
ent magnetic fields applied in-plane parallel to the iron 
strip. At zero field, the data reveal the density of states 
of the superconductor with a well-resolved energy gap 
A « 200 /ieV. With increasing magnetic field, the gap 
features broaden, and the Zeeman splitting of the den- 
sity of states for spin up and down is observed [17|. Fig- 
urel^b) shows the subgap differential conductance of the 
same junction measured at _B = (symbols) on an en- 



larged scale. Two well-resolved peaks of equal height are 
seen at V = ±23 /iV (note the expansion of the vertical 
scale by a factor of 1500). These peaks exhibit two of the 
salient features of resonant Andreev reflection via bound 
states: First, their position is symmetric with respect to 
the chemical potential of the superconductor. Second, 
the peak height is symmetric with respect to bias, as ex- 
pected for resonant Andreev reflection. 

To fit our data, we assume three contributions to the 
conductance 



9{V) ^ gTiV) + gABsiV) + go- 



(2) 



Here gT{V) describes the dominating contribution of 
« 1 mS for eV > A due to ordinary quasiparticle tunnel- 
ing, gABs{V) represents the much smaller subgap contri- 
butions due to an Andreev bound state, and 50 ~ 1 /^S 
accounts for additional subgap leakage through the thin 
oxide barrier. 

For gT(V) we have used the standard model of spin- 
polarized quasiparticle tunneling in high magnetic fields 

MM 



9t{V) 



^5^(1 TP) / n±(e)/'&. 



(3) 



where P is the spin polarization of the tunnel conduc- 
tance, /' — ~df{e + eV)/deV is the derivative of the 
Fermi function, e is the elementary charge, and the nor- 
malized density of states per spin in the superconductor 
is given by 



n±(e) 




(4) 



The complex quantities u± have to be determined from 
the implicit equation 



u± 



1 



(5) 
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where /ib is the Bohr magneton, F is the magnetic pair- 
breaking parameter and 6so = ?i/3TsoA measures the 
spin-orbit scattering strength. 

Since the observed subgap peaks were always much 
smaller than the conductance quantum, we have assumed 
a single conductance channel for (?abs(^)- The subgap 
conductance due to an Andreev bound state is then given 
by 



:fde, (6) 



2p cos (2(5 ± 61s 



where cos (5 = e/A, t'^ — t-^-t^, = r^r^, and t± and 
r± = 1 — t± are the spin-dependent transmission and 
reflection probabilities of the interface. 
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For fitting the data, we first used only the dominating 
contribution gT{V) to obtain the pair potential A, the 
magnetic pair-breaking parameter F, the normal-state 
conductance Gn, and the degree of spin polarisation P, 
following the same procedure as described in Ref. [l6| . 
To obtain a good fit of the onset of the gap features, we 
also had to adjust the temperature (T w 50 — 100 mK), 
and add a life-time broadening parameter Fit by replac- 
ing e by e + iTn in ([5]) in order to optimize the fit. For 
the spin polarization of the tunnel conductance we ob- 
tain P w 10 — 15%, typical for junctions with ultra-thin 
amorphous aluminum oxide barriers [21[. The same P 
was also extracted from non-local spin-valve experiments 
performed aX T = 4.2 K. From these experiments, we 
also obtain the coercive field of our Fe wires of about 
50 mT. 

After fixing g^iV), we then added 5abs(^) and go to 
fit the subgap peaks. We obtained typical values « 0.1 
and i_ « 10^'^ from our fits for all the observed bound 
states (note that i+ and t- are interchangeable in ([6|), 
and we have arbitrarily chosen tj^ > <_). The product 
t+i_ essentially fixes the spectral weight, while the larger 
of the two fixes the width of the conductance peaks. 
Whether t+ ^ t_ signifies a locally enhanced spin po- 
larization or simply an additional life-time broadening 
of the resonances can not be determined from the fits. 
The peak position, and consequently 16*8 1, varied strongly 
from junction to junction (see below). The result of such 
a fit, with 16*3! = 0.947r, is shown in Fig. [2jb). Upon in- 
creasing temperature, the peaks simply broaden further, 
until they disappear in the smearing of the gap features 
at about 100-150 mK (not shown). 

The density of states induced by Andreev bound states 
is spin polarized, i.e., the two conductance peaks corre- 
sponds to opposite spin. This expectation can be checked 
by observing the evolution of the peaks with magnetic 
field, as shown in Fig. [3][a). The peak features can be 
traced up to -B w 1 T, where they merge with the onset 
of the gap features given by QTiy). The magnetic-field 
dependence of the pair potential A due to orbital pair 
breaking obtained from the fits to gT{V) is negligible in 
this field range, so that the shift of the subgap peaks is 
entirely due to the Zeeman effect. As can be seen, the 
peaks shift linearly into opposite directions with a slope 
of zb/xei?, as indicated by the dotted lines. From the 
direction of the Zeeman shift, we can unambiguously as- 
sign the spin of the conductance features, as indicated 
by arrows. We can then use the spin assignment to infer 
the sign of the phase shift 9s = — 0.947r. To illustrate 
this, we have also plotted the spin-dependent density of 
states [l3l in Fig. [Stc) for the same parameters as used 
for the fit shown in Fig. [2jb). In Fig. [SUb), we display 
the data for a second contact for comparison. As can be 
seen, the energy of the bound state differs, and also the 
Zeeman shift is in the opposite direction. Consequently, 
we can assign a positive 9s = -l-0.737r. Since the coercive 
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FIG. 3. (Color online) Density of states and Zeeman shift of 
the bound states, (a) and (b) Differential subgap conductance 
of two contacts of the same sample as a function of applied 
magnetic field B. The dotted lines indicate the Zeeman shift 
of the bound states (shown for positive fields only), arrows 
indicate spin, (c) Spin-resolved density of states due to the 
bound state shown in panel (a), plotted for the parameters 
obtained from fits to (|6]), and with the sign of the spin-mixing 
angle inferred from the Zeeman shift (see text), (d) Data for 
a third contact at _B = (circles, left conductance scale) and 
in the normal state at B = 2 T (diamonds, right conductance 
scale). Both scales span 15 ^S. The lines are fits to ([2l) 
and the standard model of dynamic Coulomb blockade [22| ]. 
respectively (see text). 



field of the iron strips of about 50 mT is smaller than the 
magnetic-field increments in Fig. |3l the magnetization 
and magnetic field are always parallel for both S > and 
-B < 0. Consequently, the Zeeman shift seen for i? > 
is mirrored for i3 < 0, since spin up and down simply 
exchange their roles in the upper and lower half-plane. 
This is consistent with the view that the observed peaks 
are induced by the leakage of the exchange field of the 
ferromagnet into the superconductor. We also note that 
the observed magnetic-field dependence of the peaks ex- 
cludes the possibility that they are minigap features due 
to the proximity effect, since these should quickly disap- 
pear at higher fields. To check whether the peaks might 
be due to some inelastic processes unrelated to supercon- 
ductivity. Fig. |3Jd) compares the peaks seen in a third 
contact in the superconducting state at -B = to the nor- 
mal state data at -B = 2 T just above the critical field. In 
the normal state, a Coulomb dip is observed, which can 



be fit with the standard model [22| . The arrows indicate 
the extrapolation of the peak positions from i? = to 
_B = 2 T, using the Zeeman shift observed in the super- 
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conducting state. No peaks are seen at these positions in 
the normal state within the resolution. 

A total of 30 junctions were measured, with results 
summarized in Fig. 31 In Fig. IDJa), we show the normal- 
ized energy e+ /A of the spin- up feature, as determined 
from the fits and the Zeeman shift. In 19 of the junc- 
tions usually a single bound state was observed, and in 
some cases, there were two or three states in a single 
junction (e.g., contacts 21 and 22). 11 junctions did not 
show any subgap features, which means that there is ei- 
ther no bound state, or its energy is too close to A to be 
distinguished from the onset of the gap features. The en- 
ergy range where observation is impeded is greyed out in 
Fig. m^a). The corresponding histogram of bound-state 
energies is shown in Fig. Hljb). As can be seen, there is 
a broad distribution, with maybe a slight preference for 
values around e+/A —0.2 and 0.4, corresponding to 
6*3 « -0.97r and -|-0.757r. 

From these findings the following picture emerges: The 
simultaneous presence of strong gap features and small 
overall subgap conductance indicates that over most of 
the area of a given junction both Og and the trans- 
mission probability are small. Only for a small part 
of the junction, where the oxide barrier is thinnest, 0s 
is large, and the transmission probability is sufficiently 
large {t ^ 10"'^) to sustain Andreev transport. This 
is corroborated by the thickness variation of the oxide 
barrier seen in the HRTEM image in Fig. [Ijb). Recent 
theoretical predictions suggest that 9s can be large for an 
ultra-thin tunnel barrier with a smooth potential profile 
[l^ . The broad distribution of 0s reflects the extreme 
sensitivity of spin mixing to interface properties. 

A possible alternative scenario for the appearance of 
subgap Andreev resonances would be the presence of dis- 
crete impurity states in the disordered oxide barrier, ef- 
fectively acting as quantum dots (see, e.g., [l^llJl, and 
references therein) . This might easily explain why we see 
discrete resonances despite the fact that we have large- 
area disordered contacts. However, a couple of observa- 
tions are incompatible with this view. First, such fea- 
tures were not observed in p revious experiments on non- 
magnetic junctions [3, fabricated in the same way 
and with similar resistance-area product, where the same 
type and density of impurity states should be present. 
Second, in this case a Zeeman splitting rather than a 
simple shift of the resonances would be expected in the 
magnetic field . Third, both the energy and the width 
of the resonances would depend on the coupling strength 
to the superconducting electrode, with broader features 



at lower energy 2J]. In contrast, we did not observe 
a systematic dependence of the peak width on energy 
within the scatter of the data. And finally, magnetic im- 
purity states at the surface of a superconductor usually 
exhibit an asymmetry in peak height [2^, which we did 
not observe. Since the features are very small, one might 
also imagine a weak remnant of inelastic single-particle 
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FIG. 4. (Color online) (a) Normalized bound-state energy 
£+/ A for all measured junctions. Sample numbers are given 
at the bottom, (b) Histogram of the data from panel (a), bin 
size 0.1. 



tunneling via magnetic impurity states [26j. This would 
however be subject to the same thermal suppression as 
the overall subgap conductance, and should therefore be 
clearly visible in the normal-state conductance at high 
magnetic fields, where no features were observed at the 
expected positions in the experiment. 

In conclusion, we have shown evidence for subgap 
Andreev bound states at spin-active interfaces between 
superconductors and ferromagnets, and directly deter- 
mined the spin-dependent interface phase shift 9s, thus 
providing microscopic insight into the generation of 
triplet superconductivity at spin-active interfaces. We 
find that 6*5 can be surprisingly large, of the order of 
TT, in agreement with recent theoretical predictions for 
ultra-thin tunnel barriers. 
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